Biomes in North and South America
BiomeMap <- qtm(biomes_shp,
fill="biomes",
fill.style="fixed",
fill.labels=biomes_shp$biomes,
fill.palette=cols7,
fill.title="Biomes",
layout.legend.position=c("left","bottom"),
layout.legend.width=1.5,
layout.frame=FALSE)
BiomeMap
Figure 1. Map of North and South American biomes.
Beta Diversity Plots by Biome
BiomeBetaWeightBV <- ggplot(BiomeBetaCellsCentCov, aes(x=Type, y=Beta, fill=Type, color=Type)) +
geom_boxplot(show.legend = FALSE, fill=biome_cols_11, color="black") +
guides(x = guide_axis(angle=30)) +
#theme_minimal() + #un-comment whichever theme you want
theme_gray() +
#theme_light() +
#theme_bw() +
geom_violin(scale="count", show.legend=FALSE, fill="gray", alpha=0.35,
color="gray25") +
xlab("Biome") +
ylab("Beta Diversity") +
theme(axis.title.y = element_text(size=22),
axis.title.x = element_text(size=22),
axis.text.y = element_text(size=16),
axis.text.x = element_text(angle = 30, hjust = 1, size = 16))
BiomeBetaWeightBV

Range Size Plots by Biome
BiomeRangeBV <- ggplot(BiomeRange, aes(x=Type, y=Range, fill=Type, color=Type)) +
geom_boxplot(show.legend = FALSE, fill=biome_cols_11, color="black") +
guides(x = guide_axis(angle=30)) +
#un-comment whichever theme you want
#theme_minimal() +
theme_gray() +
#theme_light() +
#theme_bw() +
geom_violin(scale="count", show.legend=FALSE, fill="gray", alpha=0.35,
color="gray25") +
xlab("Biome") +
ylab("Range") +
theme(axis.title.y = element_text(size=22),
axis.title.x = element_text(size=22),
axis.text.y = element_text(size=20),
axis.text.x = element_text(angle = 30, hjust = 1, size = 16))
BiomeRangeBV

Alpha Diversity by Order, Split into Biomes
Orders with Maximum Alpha Diversity > 10
FacetOBRBelow10 <- ggplot(OBRBelow10DF, aes(x=Biome, y=Alpha, fill=Biome, color=Biome)) +
geom_boxplot(show.legend = FALSE, fill=biome_cols_166, color="black", outlier.size=0.7) +
#un-comment whichever theme you want
#theme_minimal() +
theme_gray() +
#theme_light() +
#theme_bw() +
#theme_dark() +
#theme_linedraw() +
#geom_violin(scale="count", show.legend=FALSE, fill="gray", alpha=0.35, color="gray25") +
xlab("Biome") +
ylab("Richness") +
theme(axis.title.y = element_text(size=22),
axis.title.x = element_text(size=22),
axis.text.y = element_text(size=15),
axis.text.x = element_text(angle = 90, hjust = 1, size = 12))+
facet_wrap(~Order,
#un-comment # of rows you want
ncol=6
#ncol=5
)
FacetOBRBelow10

Orders with Maximum Alpha Diversity between 10 and 25
FacetOBR10to25 <- ggplot(OBR10to25DF, aes(x=Biome, y=Alpha, fill=Biome, color=Biome)) +
geom_boxplot(show.legend = FALSE, fill=biome_cols_66, color="black", outlier.size=0.7) +
#un-comment whichever theme you want
#theme_minimal() +
theme_gray() +
#theme_light() +
#theme_bw() +
#theme_dark() +
#theme_linedraw() +
geom_violin(scale="count", show.legend=FALSE, fill="gray", alpha=0.35, color="gray25") +
xlab("Biome") +
ylab("Richness") +
theme(axis.title.y = element_text(size=22),
axis.title.x = element_text(size=22),
axis.text.y = element_text(size=15),
axis.text.x = element_text(angle = 45, hjust = 1, size = 12))+
facet_wrap(~Order,
#un-comment # of rows you want
ncol=3
#ncol=2
)
FacetOBR10to25

Orders with Maximum Alpha Diversity between 25 and 100
FacetOBR25to100 <- ggplot(OBR25to100DF, aes(x=Biome, y=Alpha, fill=Biome, color=Biome)) +
geom_boxplot(show.legend = FALSE, fill=biome_cols_87, color="black", outlier.size=0.7) +
#un-comment whichever theme you want
#theme_minimal() +
theme_gray() +
#theme_light() +
#theme_bw() +
#theme_dark() +
#theme_linedraw() +
geom_violin(scale="count", show.legend=FALSE, fill="gray", alpha=0.35, color="gray25") +
xlab("Biome") +
ylab("Richness") +
theme(axis.title.y = element_text(size=22),
axis.title.x = element_text(size=22),axis.text.y = element_text(size=15),
axis.text.x = element_text(angle = 30, hjust = 1, size = 10)) +
facet_wrap(~Order,
#un-comment # of rows you want
ncol=4
#ncol=2
)
FacetOBR25to100

Orders with Maximum Alpha Diversity > 100
FacetOBRAbove100 <- ggplot(OBRAbove100DF, aes(x=Biome, y=Alpha, fill=Biome, color=Biome)) +
geom_boxplot(show.legend = FALSE, fill=biome_cols_22, color="black", outlier.size=1) +
#un-comment whichever theme you want
#theme_minimal() +
theme_gray() +
#theme_light() +
#theme_bw() +
#theme_dark() +
#theme_linedraw() +
geom_violin(scale="count", show.legend=FALSE, fill="gray", alpha=0.35, color="gray25") +
xlab("Biome") +
ylab("Richness") +
theme(axis.title.y = element_text(size=22),
axis.title.x = element_text(size=22),
axis.text.y = element_text(size=15),
axis.text.x = element_text(angle = 45, hjust = 1, size = 12)) +
facet_wrap(~Order,
#un-comment # of rows you want
ncol=2
#ncol=1
)
FacetOBRAbove100

Circle Plots
Species richness for each order, grouped by biome
All bryophyte orders
circos.clear()
circos.par(start.degree = 0)
chordDiagram(CircleMatAll, grid.col = grid.col, column.col = biome_cols_11,
directional = 1, direction.type = "arrows", link.arr.type = "big.arrow",
link.arr.length = 0.05, link.largest.ontop = T, annotationTrack = c("grid"),
preAllocateTracks = 1, big.gap = 20, small.gap = 2)
circos.trackPlotRegion(track.index = 1, panel.fun = function(x, y) {
xlim = get.cell.meta.data("xlim")
ylim = get.cell.meta.data("ylim")
sector.name = get.cell.meta.data("sector.index")
circos.text(mean(xlim), ylim[1], sector.name, facing = "clockwise", niceFacing = TRUE, adj = c(0, 0.5), cex=0.5)
}, bg.border = NA)

Only moss orders
circos.clear()
circos.par(start.degree = 0)
chordDiagram(CircleMatAllMoss, grid.col = grid.col, column.col = biome_cols_11,
directional = 1, direction.type = "arrows", link.arr.type = "big.arrow",
link.arr.length = 0.05, link.largest.ontop = T, annotationTrack = c("grid"),
preAllocateTracks = 1, big.gap = 20, small.gap = 2)
circos.trackPlotRegion(track.index = 1, panel.fun = function(x, y) {
xlim = get.cell.meta.data("xlim")
ylim = get.cell.meta.data("ylim")
sector.name = get.cell.meta.data("sector.index")
circos.text(mean(xlim), ylim[1], sector.name, facing = "clockwise", niceFacing = TRUE, adj = c(0, 0.5), cex=0.5)
}, bg.border = NA)
